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Abstract 

Electromagnetic wavelets are constructed using scalar wavelets as su- 
perpotentials, together with an appropriate polarization. It is shown 
that oblate spheroidal antennas, which are ideal for their production 
and reception, can be made by deforming and merging two branch cuts. 
This determines a unique field on the interior of the spheroid which 
gives the boundary conditions for the surface charge- current density 
necessary to radiate the wavelets. These sources are computed, in- 
cluding the impulse response of the antenna. 



1 Complex distance and its branch cuts 

Electromagnetic wavelets were introduced in [K94] as localized solutions of 
Maxwell's equations. They are 'wavelets' in the historical sense of Huygens 
as well as the modern one: being generated from a single 'mother wavelet' 
by conformal transformations including translations and scaling, they form 
frames that provide analysis-synthesis schemes for general electromagnetic 
waves. Together with their scalar (acoustic) counterparts, they have been 
called physical wavelets. It was pointed out that they can, in principle, be 
emitted and absorbed causally, and applications to radar and communica- 
tions have been proposed [K96, K97, Kl] based on their remarkable ability to 
focus sharply and without sidelobes. Similar objects, long studied in the en- 
gineering literature under the name complex- source pulsed beams, have been 
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used to build beam summation methods and analyze the behavior of solu- 
tions (see [HFl] for a recent review). However, to implement the proposed 
applications to radar and communications, the wavelets must be realized 
by simulating their sources. This has proved to be difficult, and detailed 
investigations have only been made recently [HLKO, K3]. Here I present a 
new and rather complete analysis of the sources based on an insight which I 
believe is a key to their realization: the sourees can be constructed from the 
very branch cuts that give the wavelets their remarkable properties. 

Physical wavelets are based on the idea of displacing a point source to com- 
plex coordinates. Since a real translation gives nothing new, it suffices to 
discuss a point source with purely imaginary coordinates ia. It will be seen 
that this results in a real, coherent, extended source distribution parameter- 
ized by the single vector a, much as an antenna dish can be described by a 
single vector giving the orientation and radius of the dish. 

Recall the definition of the complex distance a from the imaginary source 
point ia to the real field point r, 

a{r — ia) = s/ {r — ia) ■ (r — ia) = y^r^ — — 2ia ■ r. (1) 

For each fixed source location ia ^ 0, its branch points form a circle: 

(7 = ^ r e C = {r e : r = a, a ■ r = 0}. 

It is important to note the topological difference between a ^ 0, when R^ — C 
is multiply connected, and a = 0, when C contracts to the origin and — C 
becomes simply connected. Writing 

a = p-iq, (2) 

note that (1) implies 

2 2 2 2 

r —a = p — q , a ■ r = pq, 

from which one easily obtains the relations to the cylindrical coordinates 
with z-SLxis parallel to a: 

az = pq, ap = a\/ — z"^ = \J -\- o?-\J cP- — q^. (3) 

This gives an important bound on the 'complexness' of a in terms of the 
'complexness' of its argument: 

|g| < a, or I Im (7(z)| < I Im z|, z = r — ia. 
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It follows from (3) that 



p2 

— 1 



-± ^ = 1, 

_|_ Q,2 _ g2 g2 ' 



hence the level surfaces of form a family of oblate spheroids 



Sp = {r : = const > 0} = < r : 



+ 



q2 _|_ p2 p2 



and those of form the orthogonal family of one-sheeted hyperboloids 



Tiq = {r : < q"^ = const < a^} = Ir 



a? — q^ q^ 



1 



(4) 



(5) 



To complete the picture, we include the following degenerate members of 
these families, 



So = {r 
Ho = {r 
na = {r 



< p < a, z = 0} 
a < p < CO, 2 = 0} 
/9 = 0, — oo < z < oo}. 



(6) 
(7) 



where Sq and "Hq each form a twofold cover of the indicated sets. 




Figure 1: The level surfaces of a form an oblate spheroidal coordinate system 



The families Sp and 7iq are depicted in Figure 1, along with the azimuthal 
half-plane (/)=constant. They are confocal, with the branch circle C as their 
common focal set. Note that the intersection of Sp with Tig consists of two 
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circles whose further intersection with the azimuthal half-plane consist of 
two points for each choice of {p^ , q"^ , (j)) . When p = q = 0, the two circles 
merge with the branch circle C. The set of numbers 

(cr, (f)) = (p, q, (j)), — oo < p < oo, —a < q < a, < (f) < 27t 

therefore gives a twofold cover of — C. To obtain a coordinate system, 
we must choose between the two covers, and this amounts to choosing a 
branch cut that makes a single- valued, as explained below. This will result 
in a one-to-one correspondence between (cr, (j)) and points r G not on the 
branch cut, giving an oblate spheroidal coordinate system. 

If we continue a analytically around a closed loop threading the branch circle 
C, it returns with its sign reversed. To make it single-valued, it is therefore 
necessary to choose a branch cut that prevents the completion of the loop. 
Note from (1) that 

r S> a =^ a ^ ±{r — iacosO). (8) 

The spatial region with r ^ a will be called the far zone (we need a > here 
to set the scale). Since we want a to generalize the usual positive distance 
r, we insist that 

r»o ^ p ^ r, q^acosO. (9) 

If follows from (9) that the branch cut B is bounded since it must be entirely 
contained inside any spheroid Sp with p'^ sufficiently large, and its boundary 
must be the branch circle: 

dB = C. (10) 

(The alternative is a branch cut extending from C to infinity, but this violates 
(9).) B is therefore a membrane spanning C, and any such membrane will do. 
The situation is best understood topologically. The analytic continuation of 
the distance has opened up a window connecting the two branches of r = 
^/r ■ r, thus making — C multiply connected. The spherical coordinates 
r and 9 merge analytically into cr, which is double- valued, and the choice of 
a branch cut B makes — ;B simply connected and a single-valued. 

Let (To denote the complex distance with the flat disk Sq as branch cut: 

o"o = Po - i<lo, Po > 0, -a<qo< a. (11) 



4 



The complex distance erg with B as branch cut is now defined as follows. 
Choose an arbitrary reference point ro in the far zone, where (Tb = ctq. 
To find (jg at any other point r, continue analytically along an arbitrary 
path from ro to r, with the following rule: whenever the path crosses B, 
changes sign. This gives a unique definition of ajs, and both and have 
a jump discontinuity across the interior of B. Of course, erg = ctq = on C. 

A general branch cut can be specified by a function xip, 4>) ^ 

B = {r : p = xiq, </>), -a < q < a, < ^ < 27r}, (12) 

where the cut function x must satisfy 

x{-qA) = -xiqA), x(«,27r) =x(9,0). (13) 

The first condition ensures that a changes sign across B, while the second 
ensures that it is continuous across the half-plane (f> = on each side of 
B. Note that B need not be cylindrically symmetric. We will be especially 
interested in the cuts Ba defined by the cylindrical cut functions 

Xa{q) = aSgn{q). (14) 

Note that on B^ we have 

az = pq = aq Sgn {q) = a\q\. 

If a > 0, then the values q ^ generate the upper spheroid: 

S+ = {reSa:z> 0}. 

But this does not include the branch circle C, so the bounding condition 
(10) is not satisfied. The problem is that Sgn (q) is undefined at q = 0, and 
the set of all points with g = is the degenerate hyperboloid TCq (7). Hence 
we define the part of Ba with g = as the apron bridging the gap between 
C and S^, 

B^ = {r :a<p< ^/a"^ + a"^, z = 0}. 
Thus, for a > 0, Xaiq) defines the upper spheroidal branch cut: 

Ba = 5 + U B° , dBa = C. (15) 
Similarly, x-aiq) defines the lower spheroidal cut: 

B.a=SaUB^, dB.a=C. (16) 
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As a — > 0, both cuts contract to the doubly covered flat disk spanning C, 
which is the degenerate spheroid Sq (6). 

Every branch cut B is doubly covered. Consider any simply connected, 
closed surface S containing C in its interior. Think of 5 as a balloon and 
of C as a rigid wire ring. Now deflate the balloon, and you have a branch 
cut bounded by C. In particular, if we take S = Sa and keep the upper 
spheroid rigid while deflating, the balloon stretches around the ring to cover 
the underside of and we obtain the cut (15) (see Figure 2). 

The image of a branch cut as a balloon enclosing the singular ring is similar 
to Penrose's idea of cosmic censorship in general relativity [W99, Chapter 
5], where a horizon (S) prevents the outside observer from seeing a naked 
singularity (C). In light of the connection with Newman's analytic Coulomb 
field (see the discussion below (54)), the two may in fact be closely related. 

The discontinuity of Xa in (14) causes two problems: the q = contribution 
is undefined (hence the aprons had to be chosen 'by hand'), and the 
resulting cut had a sharp edge. For computational purposes, it may be 
better to use smooth cut functions to avoid both problems. Let e > and 
define 

XM = i Im In , xiq) = aX.iq)- (17) 

vr \^ — I'q J 

For e ^ a, Xf,{q) is a smoothed version of Sgn [q) and the resulting branch 
cut closely approximates Ba without the need to define the apron separately. 
This is shown in Figure 2. 




Figure 2: The cut function x{q) = -^eiq) with e = .005 and its branch cut 
with a = 1. The two sheets have been purposely separated to show the 
double cover. 
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2 Scalar wavelets 



For any fixed choice of branch cut B, we now denote the complex distance 
simply by a. Scalar wavelets are then defined as the retarded solutions 

^(z,r) = ^^^""^^ z = r-ia, T = t-ib, (18) 

cr cr 

where we have set the propagation speed c = 1 (otherwise Qr = g{T — a / c)) 
and g is the analytic- signal transform of a driving signal go{t), defined^ as 
the convolution of with the Cauchy kernel: 



1 /■°' 



9o{t')dt' 

— , T = t-ib 

T — t' 

go{t')dt' , i {f -t)go{t')dt' 



b_ go{t')dt' i r 

27, J _^ (f/_f)2+fe2 + 27r7-< 



it' - tf + 62 

= gi{t,b)+ig2(t,b). (19) 

gi{t, b) and g2{t, b) are smoothed versions of 550 (*) and its Hilbert transform, 
with b as the smoothing parameter. We assume that go{t) decays at infinity, 
from which it follows that ^(t) is analytic in the upper and lower complex 
time half-planes C'''. The original driving signal can be recovered as the 
boundary value 

go{t) = lim [g{t - ib) - g{t + ib)] . 
0— »+o 

(The limit of the sum gives the Hilbert transform.) In particular, if ^0 
vanishes on any open interval /, this interval becomes a window between 
the upper and lower half-planes through which the functions g{t zt ib) can 
be connected so that they are both part of a single analytic function. (This 
is a special case of the edge of the wedge theorem in higher dimensions; see 
[K3].) Since every practical driving signal vanishes at least in the remote 
past, this property will be assumed. Note that this excludes time-harmonic 
driving signals, which are however idealizations. 

Now consider the numerator of (18), 

g{T - a) = g{t - p - i{b - q)). 

Suppose that \b\ < a. Then g{T — a) is undefined along the semi-hyperboloid 
where q{r) = b, except when t — p{r) is in the zero-set of go. On the other 

^This is a special Ccise of a multidimensional definition; see [K3]. 
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hand, if |6| > a, then b — q{r) vanishes nowhere and g{T — a) is analytic at 
all {r,t). Therefore we assume from now on that 



\b\ > a (20) 

so that g{T — a) is defined unambiguously everywhere. The imaginary source 
coordinates must therefore belong either to the future cone or to the past 
cone of space-time, 

b> a =^ {a,b) €V+ 
b< -a ^ (a, b) €V_ , 

which means that the complex 4- vector from the source point iy = i(a, b) to 
the field point x = {r,t) belongs either to the forward tube or the backward 
tube of complex space-time [SW64, K3], 

T^ = {z = x-iyeC'' -.yeV^}, x = {r,t), y = {a,b). (21) 

The source distribution of is now defined as a generalized function by 
applying the wave operator, 

S{z) = (5| - V^Mz) = □,*(z), (22) 

where indicates that the operator acts only on the real space-time vari- 
ables X of the field point. It is well known that 

□ ^^^"""^ = 4Trh{t)6{r) (23) 

for any differentiable function h, and this can be extended to ^(z). Since ^ 
is differentiable in r everywhere outside of the branch cut B, (23) suggests 
that 5 is a (Schwartz) distribution supported on B, a conclusion borne out by 
a rigorous analysis [K3] . The discontinuity of * across B gives a combination 
of simple and double layer terms of S on H [K3] . 

The frequency content of ^(r) determines that of * and should therefore be 
understood. Substituting the Fourier representation of into the definition 
(19) and reversing the order of integration gives 

. , I f°° dt' duj ^ . . 

duo ^ , , I e-''^^' dt' 
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The contour in the second integral can be closed in the lower half-plane if 
6 > and in the upper half-plane if 6 < 0, giving 

/oo J 
— go{i^)Q{ub)e-'^\ r = t - ib, (25) 
-oo 27r 

where &{ujb) is the Heaviside step function. Thus if 6 > 0, ^ contains 
only the positive- frequency components of qq, and if 5 < 0, it contains only 
the negative-frequency components. In either case, the factor e~'^'' in the 
extended Fourier kernel 

acts as a low-pass filter, substantially damping frequencies |a;| ^> b~^ and 
thus smoothing out ^(t). If the driving signal go is assumed real, then 
go{±io) are related by complex conjugation and therefore so are g{t^ih). If 
5o is complex, then gQ{-^uj) are unrelated and so are g{t ^ ih). 

Example: Let g{T) be the (n — l)-st derivative of the Cauchy kernel,^ 

,M = c„M = (ifl,r-5^ = |-^. (26) 

whose Fourier transform is 

/oo 
dt e'^'^Cnit - ib) = Sgn {b)@{ujb) a;"-^e-^^ (27) 
-oo 

Thus, while b acts to suppress high frequencies, n > 1 acts to suppress 
low frequencies and we end up with a band-pass filter whose effective center 
frequency and bandwidth are given by a Poisson distribution, 

= l A„ = ^. (28) 

The behavior of * in the far zone is governed by that of g{T — s). By (9), 

' 2'Ki'^{T-aY^ 2T:[{b- a cos 6) +i{t-r)Y 
is a pulse with angle-dependent duration 

T{e] = \b - a cos e\>\b\-a = r^in > 0, (29) 



^Thc driving signal is the singular distribution go{t) = (idt)" ^<5(t), but this can be 
approximated. 
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being shortest at ^ = if 6 > a and at ^ = tt if 6 < —a. While the pulse 
duration is independent of n, the strength of the peak depends jointly on 
the size of n and the smallness of 6 — a: 

MW = |5(r-a)U^«i^. (30) 

To get a measTirc of the diffraction angle, assume b > a for definiteness. Fix 
/? > and look for the angle 9p at which 

M{ep) = e-^M(O). 

Then 

(6-acose^r = e^(6-a)", 

which gives 

2 sin2(^^/2) = 1 - cos = (e^/" - 1)^^ ■ (31) 

Thus, 6p can be made small either by choosing 6 — a ^ a, or n S> /3. In 
either case, the right side gives 0^/2. A reasonable measure is obtained with 
(3 = 1. 



3 From scalar to vector wavelets 



It is well known that every electromagnetic field can be derived from a 
pair of real scalar potentials, the most well-known examples of which are 
the Whittaker and Debye superpotentials [N55]. In this section we use the 
scalar wavelet ^' as a complex Whittaker superpotential. Although this 
is equivalent to using a pair of real potentials, disentangling the real and 
imaginary parts leads to unecessarily complicated expressions, something 
like taking the real and imaginary parts of a complicated analytic function 
f{x + iy) in order to obtain two real harmonic functions. To see how bad it 
gets, note from (49) that the fields and currents contain terms of the type 
^k-3g{k) _ ^-j -yvith = 0, 1, 2. In the simplest case k = 2 (which will give 
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the radiation terms of the field), (19) gives 



Re 



gjr - cr) ^ gi{t-b,b- q) + ig2{t -b,b-q) 
a P — iq 

g{T - a)\ pgi (/ - 5. b - q) - qg^il - b. b - q) 



(T } + q^ 



Im f ~ _ Pg^ft - 6, 6 - g) + qgi{t -b,b-q) 

and it is clear that the real expressions quickly become unmanageable. Thus, 
although we work with complex potentials and fields, we view this as a very 
compact and efficient way of computing the real fields. In particular, our 
expressions contain nothing extraneous since their imaginary as well as real 
parts have a direct physical significance. This strategy is based on the 
analyticity of ^' outside the source region, which will indeed make harmonic 
pairs out the fields D and B, as seen below. 

With as a complex Whittaker superpotential, we define the retarded com- 
plex Hertz potential 



^TT I (33) 



where tt is a fixed complex polarization vector that can be seen [K3] to be a 
combination of (real) electric and magnetic dipole moments. The real and 
imaginary parts of Z 

Z = Zf^ + iZ-m 

are interpreted as electric and magnetic Hertz vectors [BW99, pp 84-85]. 
They generate a 4- vector potential (/U = 0, 1, 2, 3) by 

^0 = -V • Ze, A = dtZ^ + VxZm (34) 

which automatically satisfies the Lorenz condition 

dtAo + V • A = 0. (35) 

In turn, it follows from potential theory (or the Poincare lemma for differ- 
ential forms) that every 4-vector potential satisfying (35) can be written in 

the form (34), so this representation is quite general. (We can even dispense 
with the Lorenz condition by performing a gauge transformation on A^. See 
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[N55] for an excellent and thorough account of Hertz potentials and their 
enormous gauge group.) The real vector fields Pg and Pm defined by 

P = Pe + iPm = nZ = {D^)7v = Stt (36) 

are the electric and magnetic polarization densities. They are distributions 
supported spatially on the branch cut B. Since we are in Lorenz gauge, the 
charge-current density is = OA^, hence 

Jo = - V • Pe , J = dtPe + V X , (37) 

with charge conservation guaranteed by the Lorenz condition. The polar- 
ization densities thus act as 'potentials' for the charge-current density, a 
property inherited directly from (34). 

The reason why Hertz potentials will be so useful can be seen by computing 
the fields: 



B = VxA = VxVxZ,n + dtV x (38) 



and 



E = -VAo - dtA = VV ■ Ze - d^Ze - dtV X Zm 

= -DZe + V xV X Ze-dtV X Zm- (39) 

Taking into account (36) gives 

D = E + Pe = VxVxZe - dtV X Zm (40) 

which is a kind of 'harmonic conjugate' of (38), so the real fields D,B can 
be expressed compactly in the complex form [S41, pp 32-34] 



F = D + iB = VxVxZ + idtZ. (41) 



Note that outside the branch cut Pg = and D = E. The Hertz 
formalism thus automatically takes account of the polarization, so that the 
expression (40), if interpreted as a distribution, is valid even within a singular 
source region. 

Again I emphasize that F = D + iB is "real" in the sense that D and B 
are real, physical fields. Yet F, like is analytic in the source- free complex 
space-time region 

% = {{r-ia,t- ih) G : |5| > a, r ^ B]. (42) 
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More simply, because of their spheroidal symmetry, ^ and F are analytic 
functions of the two complex variables {a, r) in the region 

Ur, = {(a,r) G C2 : |6| > a, p / Xsiq)}, (43) 

where Xb{q) is the cut function for B. Thus D and B really are harmonic 
conjugates as suggested earlier. On the other hand, P and character- 
ize the singularities spoiling analyticity in the source region, including the 
branch points and branch cuts. This differs from the usual practice in the 
frequency domain, where (D, B) are the real parts of separate complex fields 
(-De ) Be), and it might appear that these two representations are in conflict 
since the real fields cannot be extracted by taking the real and imaginary 
parts of Dc + iB^. To clarify this, consider the frequency components of F, 

/oo 
dt e'^'F = D^ + iB^, 
-oo 

and note that since F is complex, its positive-and negative-frequency com- 
ponents are independent. Therefore a general monochromatic field consists 
of two terms, 

G^ono = e-'^^'F^ + e^'^*F_^, a; > 

= e~"-\D^ + iB^) + e"-\D_^ + iB_^) 
= e--*(r>, + iB^) + e'^\Dl + iB*J, 

where the reality conditions have been used on D^^ and B^^. The repre- 
sentation of a monochromatic field is therefore no different from that of a 
general field: 

G^mono ~ -^mono ~l~ ^-^mono 

with both fields real: 

-Dmono(i) = 2Re{e-^'^*D^} 
S^ono(i) = 2Re{e-^'^*S^}. 

Recall that for 6 > a and b < —a, the analytic signal g{T) contains only 
positive- and negative-frequency components. Therefore 

b> a =^ + iB^ = < 
h<-a ^ + iB^ = > 0. 
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This shows that the monochromatic components of the electromagnetic 
wavelets satisfy 



Bmono{t) = 2Re { 



so B trails Z) if 6 > a, and it leads Z) if 6 < —a. (Recall also that the 
pulse travels along ±a if ±b > a.) More generally, the wavelets are helicity 
eigenstates with helicity 1 if 6 > a and — 1 if 6 < —a. This concept applies not 
only the time-harmonic components but also to general time domain fields 
[K3a]. As already mentioned, using the analytic combinations of fields also 
has the great advantage of compactness and simplicity over the alternative 
of disentangling the real and imaginary parts. 

Before launching into the field computations, I want to prepare the way for 

computing equivalent currents on the spheroid Sa in the coming section. 
Let us construct this spheroid from the two branch cuts B±a given in (15) 
and (16). Let us now denote by a the complex distance with the disk Sq 
as branch cut. This will be used as a reference for defining the complex 
distance functions with B±a as cuts, which we denote by a±. Let V± be the 
volumes bounded by together with Sq, so that 



where the signs are related to the orientations of and Sq by a. The 
union and compliments of V± will be denoted by 



Now recall the rule for crossing a branch cut other than the reference cut 
Sq: aj3 changes sign. Thus, denoting the complex distance functions with 
respect to B±a by (T±, we have 



+ 



dv_ =So-s^ 



a ' 



V = V+UV_, Vl=R^-V±, V' = R^-V. 




(44) 



as shown in Figure 3. 

The field radiated jointly by the two branch cuts B±a is therefore 



'2F{a,T) inV 
Fa{a,T) = <F{-a,T)+F{a,T) mV+ 
^F{a,T) + F{-a,T) mV_. 
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Figure 3: Values of the branches cj± of the complex distance function deter- 
mined by the branch cuts B±a, given in terms of the branch a determined 
by the disk Sq. 



Observe that there is no field discontinuity in going from to V_ , hence 

FJa,T) = < ^ ' ' (45) 
^ ^ \F{-a,T)+F{a,T) inV ^ ^ 

as depicted in Figure 4. 

The transition a — ^ —a across a branch cut turns retarded fields into ad- 
vanced fields since 

*(-a,r) = -^^^. (46) 
a 

Although advanced fields are usually associated with acausal behavior, there 
is a perfectly causal explanation for (46). Consider the field radiated back- 
ward from Ba, as observed in V_f_. Due to the curvature of the back side of Ba, 
this field converges toward the the focal ring C and, having passed though, 
it is no longer in and therefore diverges normally. A similar argument 
explains why the field emitted forward from B-a first converges toward C 
and then diverges away from it. The usual acausal behavior associated with 
advanced fields is due the the assumption that they remain advanced for 
the indefinite future. (This argument also applies to time-reversed acoustics 
[FO], where time reversal occurs only in a bounded space-time region.) 
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Figure 4: Interior and exterior fields radiated by the oblate spheroid Sa, 
represented as a combination of the two branch cuts 13-i-a. 

It was shown in [K3] that the sources of ^'(cr, r) and cr, r) are equal and 
opposite; that is, they form a source-sink pair: 

D-^{-a, t) = -□*((7, r) = -S. (47) 

The proof is trivial for real point sources, where 

□ ^^^^^ = -4Trgo{t ± r)S{r) = -4Trgoit)6{r). 

But it is more subtle for complex point sources because the extended delta 
function 

S{z) = — , z = r — ia 

with a 7^ fixed, is not supported at a single point but on the entire branch 
cut B and thererfore 

f{a)~6{z) + /(0)5~(z). 

In fact, the left side is not even defined since a is discontinuous precisely on 
the disk supporting 6{z); therefore, some care must be used in proving (47). 
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Equation (47) shows that the interior superpotential *((7, r) + cr, r) is 
sourceless, as are the Hertz potentials and electromagnetic fields derived from 
it. The interior field is 



Fo((7,T)=F(c7,r)+F(-(7,T). 



(48) 



Let us first compute the exterior field which will give the interior field 
by symmetrizing with respect to a. Let 



9r 9r . o / \ 
2 ' 9r = Otgr{T - a) 



a 



u = V(J = Vp — iVq = — 

a 

and note that w is a complex unit vector: 

z ■ z 



9r 



2^' 



a 



u - u 



I. 



a'- 



Thus 



V X Z = V* X TT = X TT 

V X V X Z = X (u X tt) + *'V X (w X tt), 
and by a simple computation, 

W X (w X tt) = — TT, A = tl • TT 
Alt + TT 



V X (w X tt) = — 



Therefore (41) gives 



F = V(Aw - tt) (Aw + tt) + i^'u x tt 

a 



or 



^ + % + %lAn- 



9r 9r 9r 



TT — Z 



9r ^gr 



This expression will be written compactly as 



L\u — Mtt — iNu X TT 



W X TT. 

(49) 
(50) 
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where 



L = 


9r 


+ 


3gr 




a 




(72 


M = 


9r 


+ 


9r 




a 




a2 ■ 


N = 


9r 


+ 


gr 




a 




(72 • 



^ (73 
(73 

(51) 



We now examine the far field to see under what conditions the polarization 
vector TT gives the strongest beams. In the far zone (9) we have 



Therefore 



r S> o ^ a^r — iacosO, tt « e^. 



Ffar = ( Ae^ — 7T — tCr X 7T) 

r 



where 



— ^(TTx + ie^ X TTx), (52) 



7r_L = TT — (tt • er)er 



is the component of tt orthogonal to r which, as expected, is the only one 
that matters in the far zone. Note that while we have replaced (7 by r in 
the denominator of (52), the presence of Im a ~ — rcos^ in g{T — a) plays 
an essential role in determining both the collimation of the beam and the 
duration of the pulse, as already seen in (30) for ^(r) = C„(r). 

The far field satisfies the helicity condition 

iSr X Ffar = Ffar , (53) 

or equivalently 

-^far = S-r X Dfai, -Dfar = — 6^ X Sfar . 

As we are interested mainly in the paraxial region of the far zone, the most 
efficient choice of tt is orthogonal to a. Since 

r — ia p z — ia ^ 
r = p + za u = = — I a, 

(7 (7 (7 
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this implies 



A = tt • TT = 



P ■ TT 



The far-zone energy density is 




and, by (53), the far-zone Poynting vector is 




:ar 




since Ffgj- is orthogonal to e^. 

4 Equivalent currents 

In principle, the scalar source generates the charge-current density by (36) 
and (37). But this would involve not only the messy disentangling of the real 
and imaginary parts of P = Sir (with both factors complex), but also deal- 
ing with the singular nature of S. While S is well-defined mathematically 
as a distribution [K3], it seems to be of little direct value from a practical 
point of view. Since S is supported on the branch cut B, one expects the 
electromagnetic sources to consist of a surface charge density jo and a sur- 
face current density j. But it turns out that these these surface sources are 
singular on the branch circle C, where a = 0. The essence of the problem 
can be understood from a careful analysis, given in [Kla], of a much simpler 
case, which we now recall. 

Example: The analytically continued Coulomb field due to point charge of 
strength Q = 1 is 



Newman [N73] has shown that this can be identified with a real electromag- 
netic field {D, B) by 



C{r - ia) = -V 



1 



r — 



(54) 



C = D + iB 



(55) 
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interpreted as the flat-spacetime (zero-mass) limit of the Maxwell field in 

the Kerr-Newman solution in general relativity [N65], which represents a 
spinning black hole of unit charge.^ It is instructive to compute the surface 
sources on a branch cut, which for simplicity we now take to be the disk Sq 
defined in (6) . On the upper and lower faces of <So we have 



p —>■ +0, z —>■ ±0, a cfi — p^, 

hence 

fp a 



^47r(a2 - p2)3/2 ^ 47r(a2 _ ^2)3/2 

The jumps in D and B across the cut are therefore 
6D = ..A^, SB 



27r(a2 - p2)3/2 ' 27r(a2 - p2)3/2 • 

Since 5D is orthogonal and 5B is tangent to 5o, it follows that the magnetic 
surface charge- and current densities vanish as required. The electric surface 
densities are given by [J99, p 18] 

j, = -a.SD = -^-^^^^ (56) 

^ „ cax p cpes 

J = a X SB — 



27r(a2 - p2)3/2 27r(a2 - p2)3/2 ' 

where we have inserted the speed of light (taken ealier to be c = 1) for 
dimensional reasons. Before discussing the problem with (56), note that if 
we define the local charge velocity by 

v = i = 2Pe^, (57) 
a a 

its linearity in p suggests a 'hydrodynamic' interpretation of Sq as a rigidly 
spinning charged disk with angular velocity 

n = c/a. (58) 

In particular, the rim C is moving at the speed of light. While this con- 
clusion seems bizarre in ordinary electrodynamics, it is entirely consistent 

^ When D in (55) is reinterpreted as a Newtonian force field, then S is a gravitomagnetic 
field related to the 'dragging' of Einsteinian spacetime in the vicinity of a spinning body. 
Evidently, this effect survives the flat-spacetime limit as the conjugate-harmonic partner 
to Newtonian gravitation. 
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with the origin of the field C as the residual Maxwell field of a charged, 
spinning black hole. The investigation in [Kla] has sparked a renewed in- 
terest in Newman's original paper [N73], leading to similar interpretations 
of linearized gravitational fields [N2] and a generzlized Lienard-Wiechert 
field where the radiating point source moves along an arbitrary trajectory 
in complex spacetime [N4]. Our antennas will be similar, but their source 
is a dipole following a complex trajectory and not a monopole, and so their 
charge-current densities arc generated by polarizations. 

We now come to the main lesson taught by this example. C is an analytic 
continuation of the Coulomb field of a point source with charge Q = 1. If 
the continuation is to make physical sense, the total charge should remain 
unchanged. This is contradicted by j'q, which is not only strictly negative 
but whose total charge on Sq is — oo! To resolve this difficulty, it is necessary 
to treat the charge-current density as a singular volume distribution, just as 
the scalar source S = was treated in [K3] . The inhomogeneous Maxwell 
equations now state that the (volume) charge- and current density are 



while the homogeneous Maxwell equations require that Jq and J be real. 
Taken as definitions of the sources in the sense of generalized functions, it 
was shown in [K3a] that (59) indeed give a sensible answer. The equivalent 
surface sources on any spheroid Sa with a > are defined by 



where the outgoing unit normal Sp on Sa is computed in the Appendix. 
These sources are found to be com,plex, which means that they include a 
magnetic charge-current density; the latter vanishes in the limit a — > 0, in 
agreement with the above conclusion. The advantage of using a > is that 
the sources are smooth and hounded, with a total charge Q = 1 as required. 
As a ^ 0, they decompose into surface sources on the interior of the disk 
which coincide with (56), plus line sources on the rim C. The line sources 
carry a total charge of oo, but when the entire source distribution is treated 
as a generalized function, it carries the correct total charge Q = I. The 
problem with (56) is that the jump conditions (using infinitesimal pillboxes 
and loops) can be applied only on the interior of the disk and not on its 
boundary C. A similar argument applies to every branch cut, showing that 
caution must be exercised in computing equivalent sources, a lesson we will 
recall when computing the currents required to produce electromagnetic 
wavelets. 



Jo = V • F, J = -dtF -iVxF, 



(59) 



jo = Bp ■ SF, j = -iep X 5F, 



(60) 
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Finally, we turn to computing the equivalent sources for F on the spheroid 

Sa- Some important properties of equivalent real scalar surface sources were 
analyzed in [HLKO], but their connection to the vector case and, specifically, 
to our topological use of branch cuts, remains to be explored. 

According to (45) and (48), the jump in the field across the spheroid is 



6F{a,T) = F{a,T)-F{-a,T) (61) 



where the complex distance a with respect to Sq is continuous across Sa- 
Unlike the sum (48) of retarded and advanced fields, the difference (61) 
does have sources and they are confined to the surface Sa, which we shall 
presently compute. Begin by writing (50) in the more explicit form 

F{a, t) = L(cr, t)\u - M{a, T)7r - iN{a, t)u x tt (62) 

with L, M, N given in terms of the retarded signal gricr, r) = g{T — a) by 



Lia, 


r) 


_ 9r 

a 


+ 




M{a, 


r) 


_ 9r 

a 


+ 


9r ^ 9r 

a"2 


N{a, 


r) 


_ 9r 

a 


+ 


9r 
(72 • 



Define the mixed signals g± by 

g±{<^, r) = g{T -a)± g{T + a) 

and note that 

a — — £7 w = > — tt, A = w • TT — —A. 

(7 

Then we obtain the following expression for the field discontinuity: 

6F = L{a, t)Xu - M{a, t)tv - iN{a, t)u x tt, (63) 

where 



L((7, 


r) = 


'9± 
a 


+ 


3g- ^ 3g+ 

£72 (7^ 


M((7, 


r) = 


h. 

a 


+ 


0-2 (7^ 


iV((7, 


r) = 


9_ 

(7 


+ 


9+ 
(72 • 
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Before going on to compute the currents, note that (45) can be modified so 
that the interior field is any source-free solution of Maxwell's equations, ie, 




2F{a, t) in V 
Fint(r,t) inV 



The choice of an interior solution other than Fq (a, r) will of course modify 
the equivalent sources on Sa- However, unless Fint fits into the spheroidal 
geometry, the resulting sources can be expected to be much more com- 
plicated and unnatural, and probably will not benefit from the 'magic' of 
complex source points. Probably the most general class of interior fields 
that do fit the geometry consists of arbitrary multiples of Fq, ie, 



and all computations below easily generalize to this case. However, only 
when = u = 1 can the radiating surface he interpreted as a combination of 
branch cuts! I believe that this case is the most natural and expect it also 
to be the most useful. For this reason, only it will be treated here, although 
our results easily extend to the case (64). 

I now compute or estimate the various inner and outer products needed in 
(60) . Free use will be made of the results derived in the Appendix, and there 
is no pretense of rigor. I will assume that 



F;^t{r,t) = vFo{a,T). 



Then (61) is replaced by 



(5F((7,t) =/xF((7,t) -i/F(-(7,T), IJL + U = 2, 



(64) 



< a < a. 



which means that the spheroid is rather flat. By (3) and (75), 



pK, ^a^ -(^, z = aq, |Vp| — , |Vg| « — . 



Thus 



u = 



1 — 



p+ (z — ia)a p — ia 



(65) 



a 



a 



Recall that tt is orthogonal to a, so that 
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and thus 



^ r ■ TT p ■ TT pTTo 
X = U ■ TT = = = -. 

a a a 

Prom (76) in the Appendix, the outgoing unit normal on Sa is 

Vp ar + qa qa q ^ 

Bp = —— = , - , = ~ — , = -r—; a. (do) 

^ |Vp| + gVa^ + ^2 + ^2 |a| 

The term has been retained in the denominator to control the singularity 
at the equator. (This is the main advantage of using Sa instead of So-) The 
approximation (66) fails very near the equator q = 0, where Bp is far from 
parallel to a, but the analysis in [HLKO] suggests that, for scalar wavelets at 
least, the immediate vicinity of g = can be ignored. More precisely, it was 
shown that for time-harmonic driving signals of frequency co, the effective 
aperture, emitting most of the radiation, consists of the front surface of the 
disk Sq parameterized by 

A;~^<g<a, ie, p'^ < - 1/k'^, k = uj/c. (67) 

Of course, this has significance only if fca > 1. Lower frequencies generate 
mostly a reactive field that swirls around the source region and is eventually 
reabsorbed.^ Thus, to obtain a high radiation efficiency, it is necessary to use 
signals g{T) with little low-frequency content, such as linear combinations of 
high-order derivatives of the Cauchy kernel (26). (Of course, a careful rep- 
etition of the analysis needs to be made specifically for the electromagnetic 
case.) 

The inner products needed to find jo are 

Bp - u = Bp - {Vp — iVq) = \Vp\ ~ 

|(t| 

Bp - -K ^ - — ra-7r = 

\a\ 

Bp ■ {U X 7V) = TV ■ {Bp X U) = TV ■ [Sp X (Vp — iV q)) 

= i|Vg|7r • {Bq X Bp) = i|Vg|7r - b^'^ -— tt^. 



*This is consistent with the general fact that 'DC components do not propagate.' It is 

also the basis of one of the close connections between electromagnetic wavelets and math- 
ematical wavelet theory, since it amounts to an admissibility condition on electromagnetic 



wavelets [K94, p 214]. 
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The outer products needed for j are 

epXu = ^|Vg|e^ « i-^ 
\a\ 

Cp X TT ~ j — j- a X TT ~ -j — j- [TTpefj) — TTfpP) 
Bp X (u X Tt) = (Cp ■ 7V)U — (ep ■ U)7V — (Cp • u)7r — j — j- TT. 

Using these in (60) gives the approximate surface charge density 



c|c|jo ~ LapiTp + NapTTcf) (68) 



and the approximate surface current density 



a\a\j m [tp^T^p - Ma^TTp + Napn^^ + [mg'^tx^ + Naa-K^ Bp. 



(69) 



As expected from our example of the analytic Coulomb potential, the equiv- 
alent sources on a spheroid are complex, indicating the presence of unrealiz- 
able magnetic charges. Since the magnetic sources in that example vanished 
as a — 0, it is reasonable to hope that this will also be the case here. As 
the spheroid Sa with < a ^ a is very flat, it may be possible to choose 
the phase of the polarization vector tt (representing the mixture of electric 
and magnetic dipoles) so as to minimize the magnetic sources over Sa ex- 
cluding the vicinity of the rim C, and the latter region can be ignored for 
highly oscillatory driving signals as shown in [HLKO]. This question will be 
addressed in detail elsewhere. 

Finally, we compute the im,pulse response of the antenna, ie, the sources 
when the driving signal is the impulse 

go{t)=d{t) => 5(r) = ^ = Ci(r). 

Notice that the real point source version of the scalar wavelet (18) is then 
the retarded propagator for the wave equation, 

* ^ *o {r, t) = ^^^^ , = 4TrS{r , t) , (70) 

r 

where the precise relation between * and is given in terms of complex- 
distance potential theory in [K3] . The mixed signals are 

1 1 r , 2 2 

S+ = 7r~f ^ + o- ^ _L \ = ~ — ' ^^^^^ u = T ~a 

ziTT^T — a) 2nr{T + a) ittu 

_ 1 ^ _ ^ 

^ 2m{T — a) 2m{T + o) inu 
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and their time derivatives are 



This gives 



9+ — '■ 2~' 9+ — '■ Q 

2aT .. 6rV + 2c73 
9- = --■ — 2' 9- = ■■ — 3 ■ 



^ _ 15a W - 10a- r' + 

^ = 

~ _ 3o^+6aV-T^ 

which can be substituted into (68) and (69) to obtain the impulse response. 

In view of the discussion following (67), we are actually more interested in 
the system's response to the bandbass signal in (26) , 

Cn{r) = m--'C,{r) = = {-d,T-'C,{r). 

(n) 

The induced surface source can be computed directly from the impulse 
response: 

ji-) = {-d,r-'j,. (71) 



5 Concluding note 

Source-free relativistic fields always extend analytically to the double tube 
domain T± (21) of complex space-time, as explained in [K3]. I find it quite 
remarkable that the extension ^(cj, r) of the propagator (70) generates fields 
with spatially compact sources that are analytic in the source-free parts 7^ 
of complex space-time obtained by removing the world tubes swept out by 
the sources. The boundary values of these analytic fields then characterize 
the singular sources, as shown above. 
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6 Appendix 



The complex unit vector u is given by 



It = Vcr = - = Vp - iVq, (72) 



hence 



pr + qa pa - qr 

Vp = 2 _L 2 ' "g = 2 _L 2 • '^^ 



Note that 

u . u = 1 |Vp|2 - |Vgp = 1, Vp ■ = (74) 

and 

p^ + p + 2 + 

which gives 

^2 I „2 „2 Jl 



2 r + a IV7 |2 Q -g /v.N 

= 2 I 2 ' ^g = 2_L 2 - (^^) 
p2 _|_ p2 _|_ ql 



The unit vectors in the directions of increasing p and q are therefore 

+ qa 

— / / 

^p2 _|_ q2^p2 _|_ q2 

pa — gr 



(76) 



eg 



^p2 _|_ q2^a2 _ q2 



Acknowledgements 

It is a pleasure to thank Drs. Richard Albanese, Iwo Bialynicki-Birula, 
Ehud Heyman, Ted Newman, Ivor Robinson, Andzej Trautman and Arthur 
Yaghjian for friendly discussions and suggestions related to this work, and 
David Park for generous help with the figures using his DrawGrahics pack- 
age. I am especially grateful to Dr. Arje Nachman for his sustained support 
of my research, most recently through AFOSR Grant #F49620-01-1-0271. 



27 



References 



[BW99] M Born and E Wolf, Principles of Optics, seventh edition. Cam- 
bridge University Press, 1999 

[FO] M Fink D Cassereau, A Dcrode, et al.. Time-reversed acoustics. Rep. 
Prog. Phys. 63:1933-1995, 2000 

[HFl] E Heyman and L B Felsen, Gaussian beam and pulsed beam dy- 
namics: Complex source and spectrum formulations within and beyond 
paraxial asymptotics. Journal of the Optical Society of America 18:1588- 
1611, 2001 

[HLKO] E Heyman, V Lomakin and G Kaiser, Physical source realization of 
complex-source pulsed beams, .Journal of the Acoustical Society of Amer- 
ica 107:1880-1891, 2000. http://www.wavelets.com/OJASA.pdf 

[J99] J D Jackson, Classical Electrodynamics, third edition. John Wiley &; 
Sons, New York, 1999 

[K94] G Kaiser, A Friendly Guide to Wavelets. Birkhauser, Boston, 1994 
(sixth printing, 1999) 

[K96] G Kaiser, Physical wavelets and radar: A variational approach to re- 
mote sensing, IEEE Antennas and Propagation Magazine 38 #1, pp. 15- 
24, 1996. 

http: / /www.wavelets.com/96AP.pdf 

[K97] G Kaiser, Short-pulse radar via electromagnetic wavelets, in 
f/Jtra- Wideband, Short-Pulse Electromagnetics 3, pp. 321-326. C 
E Baum, L Carin and A P Stone, eds.. Plenum Press, 1997. 
http: / /www.wavelets.com/96AMEREM.pdf 

[KO] G Kaiser, Complex-distance potential theory and hyperbolic equations, 
in Clifford Analysis, J Ryan and W Sprossig (editors) Birkhauser, Boston, 
2000. http://arxiv.org/abs/math-ph/9908031 

[Kl] G Kaiser, Communications via holomorphic Green functions. In Clif- 
ford analysis and Its Applications, F Brackx, J S R Chisholm and 
V Soucek, editors. Plenum Press, 2001. http://arxiv.org/abs/math- 
ph/0108006 



28 



[Kla] G Kaiser, Distributional Sources for Newman's Holomorphic Coulomb 
Field. Invited paper, Workshop on Canonical and Quantum Gravity III, 
Banacli Center, Polish Academy of Sciences, Warsaw, June 7-19, 2001. 
http: / / arxiv.org/ abs / gr-qc /0108041 

[K3] G Kaiser, Physical wavelets and their sources: Real physics in complex 
space-time. Topical Review, Journal of Physics A: Mathematical and Gen- 
eral\o\. 36 No. 30, R291-R338, 2003. http://www.iop.org/EJ/toc/0305- 
4470/36/30 

[K3a] G Kaiser, Helicity, polarization, and Riemann-Silberstein vortices. 
To appear in Journal of Optics A, Special Issue on Singular Optics, 
http: / / arxiv.org/ abs/math-ph/0309010 

[N65] E T Newman, E C Couch, K Chinnapared, A Exton, A Prakash, and 
R Torrence, Metric of a rotating, charged mass. Journal of Mathematical 
Physics 6:918-919, 1965 

[N73] E T Newman, Maxwell's equations and complex Minkowski space. 
Journal of Mathematical Physics 14:102-103, 1973 

[N2] E T Newman, On a classical, geometric origin of magnetic moments, 
spin-angular momentum and the Dirac gyromagnctic ratio. Physical Re- 
view D 65:104005, 2002. http://arxiv.org/abs/gr-qc/0201055 

[N4] E T Newman, Maxwell Fields and Shear-Free Null Geodesic Gongru- 
ences, preprint, January 2004 

[N55] A Nisbet, Hertzian electromagnetic potentials and associated gauge 
transformations. Proceedings of the Royal Society of London A231:250- 
263, 1955 

[S41] J A Stratton, Electromagnetic Theory. McGraw-Hill, New York, 1941. 

[SW64] R F Streater and A S Wightman, PCT, Spin and Statistics, and All 
That. Addison- Wesley, 1964 

[W99] R Wald, Black Holes and Relativistic Stars. University of Chicago 
Press, 1999 



29 



